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Yu  Wenyan 

(Shanghai  Institute  of  Optics  and  Fine  Mechanics,  Academia  Sinica) 

Abstract 


The  gain  of  a  high  power  multi-pass  amplifier  is  studied. 

The  conditions  for  raising  the  energy  extracting  efficiency  are 
discussed.  We  obtain  the  general  rule  to  select  the  optimum 
parameters.  This  result  has  been  used  to  build  a  disk  amplifier; 
there  are  three  optical  path;--  in  the  amplifier;  and  the  laser 
beam  of  100  mm  aperture  passes  through  it  six  times,  back  and 
forth.  Theoretical  calculations  are  in  agreement  with  the  experi¬ 
mental  results. 

I.  THE  ENERGY  AND  GAIN  OF  A  HIGH  POWER  MULTI-PASS  AMPLIFIER 

The  operation  of  a  laser  beam  passing  through  an  ampli¬ 
fier  several  times  is  called  multi-pass  amplification.  The  pur¬ 
pose  of  studying  multi-pass  amplification  is  threefold.  Firstly, 
the  operation  can  effectively  upgrade  the  efficiency  of  the  ampli¬ 
fier.  Secondly,  in  a  laser-induced  nuclear  fusion  experiment, 
the  target-reflected  laser  will  return  to  the  amplification  system. 
Therefore,  it  is  very  important  to  evaluate  the  gain  and  energy 
of  the  reflected  laser.  This  is  a  typical  double-pass  amplification 
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problem.  Thirdly,  the  amplification  of  a  pulse  array  can  be 
treated  essentially  as  a  multiple-pass  problem. 

The  laser  pulse  being  studied  operates  in  the  nano¬ 
second  and  sub-nanosecond  range.  The  classical  theory  and 
simplified  three  energy  level  model  can  be  used  to  calculate 
the  gain.  Then,  the  variation  of  energy  density  of  the  laser 
beam  in  the  amplifier  is  determined  by  the  following  equation: 

^~£./9u[l-<r£/lv)-«£,  (1) 

2 

where  E(Z)  is  the  energy  density  of  the  laser  beam  (joules/cm  ); 

B  (Z)  =  o  N  (Z)  is  the  initial  gain  coefficient  (cm-'L);  E  =  hvl'ia 
o  2  s 

is  the  saturated  energy  density  ( joules/cin  ) ;  a  is  the  absorption 

coefficent  (cm-1). 

This  is  a  non-linear  equation  which  must  in  general  be 
solved  with  the  help  of  a  computer.  We  have  been  observing  the 
effects  of  light  absorption,  transverse  isotropy  and  light  beam 
divergence  on  the  gain.  Based  on  those  effects,  we  determined  the 
accuracy  and  application  range  for  many  approximate  analytical 
formulas,  such  as  "concentrated  loss  approximation",  "uniform 
distribution  approximation",  "equi-ef ficiency  absorption  coeffic¬ 
ient,"  etc.  This  article  will  concentrate  on  the  discussion  of 
the  gain  problem  of  multi-pass  amplification. 

In  order  to  calculate  the  gain  of  a  light  beam  passing  the 
amplifier  each  time,  the  key  problem  will  be  to  understand  fully 
the  variation  rule  for  the  density  of  reversed  particles.  This 
is  different  from  the  case  of  single-pass  amplification.  It  is 
not  only  controlled  by  the  excited  transition,  but  it  is  also 
related  to  the  relaxation  effect  between  energy  levels.  The 
excited  transitions  dominate  at  the  moment  when  light  beams  pass 
the  amplifier,  while  during  the  time  interval  between  the  departure 


of  the  prior  laser  and  the  arrival  of  the  incoming  laser’ ,  the 
relaxation  effect  is  a  decisive  factor. 


We  first  discuss  the  sign! fi cance  of  the  relaxation 

effect.  In  neodymium  glass,  there  are  two  categories  for  the 

relaxation  between  energy  levels:  one  is  the  "thermalization" 

4  4 

of  sublevels,  metastable  state  (  F0/-)  and  final  state  (  I..,  .  ), 

3  /  <=  11  / z. 

i.e.  the  particle  spontaneously  tends  to  reach  the  thermal  equi¬ 
librium  distribution;  the  other  category  is  the  vacant  final 

4 

state,  i.e.  the  particle  in  the  state  spontaneously  trans¬ 

fers  to  the  ground  state.  The  energy  level  structures  of  three 
types  of  neodymium  glasses  manufactured  in  China  are  shown  in 
Figure  1.  In  thermal  equilibrium,  the  number  of  particles  at  a 
sublevel  follows  the  Boltzmann  distribution.  For  a  metastable 
state,  the  following  relation  is  obtained: 

£4*»0.434o  (?) 


For  the  energy  level  of  a  final  state,  the  relation  is: 


]C\ 

0.619 +0.404 +  0.286 
+0.202  +  0.140  =1.65lc 


(3) 


When  the  particle  distribution  deviates  from  the  equilib¬ 
rium  state,  under  certain  external  influences,  the  relaxation  will 

restore  the  equilibrium.  The  rate  of  the  recovery  can  be  described 

phenomenologically  by  a  characteristic  time,  i.e.  the  thermalization 
time  of  the  sublevel  T.  The  value  of  T  is  roughly  on  the  order  of 

10  ns  [1],  On  the  other  hand,  the  evacuation  rate  of  the  particle 

number  at  the  final  state  level  is  described  by  the  lifetime  of 
the  final  state  t.  Various  measurements  on  i  have  been  reported, 
generally  in  the  range  of  ten  to  one  hundred  nano-seconds  [2]. 


)!577  «-m-’ 


Figure  1.  Energy  level  structure  of  the  N 
neodymium  glass.  0312 

As  a  result,  the  particle  density  at  any  sublevel  will  vary  with 
time  in  accordance  with  the  following  equations. 

For  metastable  state: 


n(<)  =n+  [n(0)  —  n]e  t,T , 
y 


(4) 


For  final  state: 


m(<)  *»  {m  +  [m(0)  —  vi]e~,/T}e~,/r0  ( 5 ) 

where  n(0),  m(0)  are  particle  densities  at  t  -  0,  n,  m  are  the 
values  at  thermal  equilibrium,  which  are  determined  by  Equations 
(2)  and  (3).  The  lifetime  of  the  metastable  state  is  not  taken 
into  account  in  Equation  (4),  simply  because  it  lasts  for  several 
hundreds  of  micro-seconds  and  does  not  affect  the  multi-pass 
amplification  process  significantly. 

It  is  easy  to  calculate  the  variation  of  particle  distri¬ 
bution  induced  by  excited  transitions.  Assume  the  initial  reversed 
particle  density  to  be  N  ,  then  the  density  of  residual  reversed 
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particles  after  amplification  can  easily  be  derived  by  the  speed 
equation,  leading  to: 


(6) 

Thus  the  particle  density  at  upper  energy  level  has  been  reduced 
by  A : 


No-N 

2 


N0 


[1 


(7) 


Combining  the  results  of  excited  transitions  and  relaxa¬ 
tion  effect,  a  general  rule  for  variation  of  particle  density 
during  the  entire  process  of  multi-pass  amplification  can  be  set 
up.  Suppose  that  the  particle  density  distribution  in  material 
is  described  by  the  following  four  quantities  prior  to  the  k-th 
amplification: 

N  1  —  the  difference  of  particle  densities  at  upper 

and  lower  working  energy  levels,  (n’-m’); 

S'  —  the  sum  of  particle  densities  at  upper  and  lower 
working  energy  levels,  (n'+m'); 

82'  —  the  sum  of  particle  density  at  each  sub-level  in 

the  metastable  state; 

S^'  --  the  sum  of  particle  density  at  each  sub-level 

in  the  final  state. 

The  particle  densities,  after  consecutive  processes  of 
amplification  and  relaxation  and  before  the  (k+l)-th  amplifica¬ 
tion,  will  be  changed  respectively  to: 

AC  (3),  *:.(«).  *1(8),  «I(8). 
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where  5  is  the  time  interval  between  the  two  amplifications. 
According  to  the  formula  previously  derived  and  referring  to 
Figure  2,  the  following  relations  are  obtained: 


A-u'h  [/£  (t +e-*<')  +  -fa  ’)  (1  -  e"'T) 

-A'OU  +  A.e-^), 

*"(«)-  ')  +4(1+^  ]^'+[-4r+_Ar^.](1_.-^ 

-J'C  A,-Aie-»'), 

*5(8)*  it- a, 

*? (<5)  --  (s'-* 


(8) 


where 


l  +  k  ’ 


5 


^4a 


l  +  A-a  ’ 


y 


and  a'  is  given  in  Equation  (6). 
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-n* 

#T- «“(«) 

*  n* 

-a 

ri'ils  -n "*  {ri  -  &  -  S”)r*'r 

-  m*  _ 

e 

jC 

m'  _ _ 

w’  ♦  s 

»"  (5)  -  Cm"  +  (m’  ♦  £  -  m ') 

(A) 

<£> 

(O 

Figure  2.  Variation  of  particle  distributions:  A)  initial 
distribution;  B)  variation  induced  by  excited  transitions; 

C)  variation  induced  by  relaxation  after  a  time  interval. 

This  is  a  set  of  recurrence  relations.  Once  the  initial 
particle  density  and  input  laser  energy  density  at  each  pass  are 
given,  the  reversed  particle  density  and  energy  gain  after  any 
amplification  can  be  calculated.  Equation  (8)  is  a  general 
formula  which  can  be  applied  to  any  situation  in  practical  oper¬ 
ations.  The  formula  also  covers  many  special  results  derived  in 


previous  articles. 


Since  the  practical  implication  of  the  double-pass 
amplification  is  more  extensive  than  others,  and  the  approximate 
solution  can  be  obtained,  we  will  discuss  double-pass  amplifica¬ 
tion  in  more  detail. 

It  can  reasonably  be  assumed  that  the  initial  reversed 
particle  density  in  the  operating  material  is  uniformly  distributed 
before  the  first  laser  train  enters  the  amplifier.  In  addition, 
the  final  state  is  nearly  evacuated,  and  the  energy  levels  in  the 
metastable  state  are  thermalized .  Then,  the  initial  conditions 
will  be: 


*n*=Arn,  “0,  82  «•  (1  +  ij)iVo  (9) 

The  recurrence  relation.  Equation  (8),  can  be  amplified  to: 


K( d)~N'0-(A!l+A16-*'')A', 

*S(8)  -  K-  (da-yl1«“9/,)d', 

*It«)  -  A'e-S'\ 


Thus,  the  gains  in  the  first  round  and  the  second  round  can  be 
represented  by  the  following  expressions: 


O’(L) 


E'(L) 

K 


"'0 


Me¬ 


dia) 


0"(L) 


_  E"  (L) 

K 

E*-l)yo]p-al0 


(11b) 
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where 

second 


E 


-  ,  E 

o  '  o 

passes 


light  at  the 
signal  gains 
defined  as: 


"  are  energy  densities  of  input  light  at  first  and 
,  respectively;  E'(L),  E"(L)  are  energy  densities  of 
output  terminal  of  the  amplifier;  go ' ,  go"  are  small 
in  the  first  and  second  amplifications,  which  are 


“9o^p\~-[E'(L)~E’0 

+«\h0E\Z)dZ]}o 


(12a) 


(12b) 


Equations  (10),  (11),  (12)  are  complete  solutions  for  the 
double-pass  amplification  problem.  The  parameter  A=A^  +  A1e~>/r 
represents  the  influence  of  the  relaxation  effect,  which  varies 
with  the  time  interval  §  .  The  following  table  lists  several 

typical  cases. 


« 

n 

T<3<  x 

Relaxation 

degree 

Thermalizat ion 
and  evacuation 
are  neglected 

Basic  therma- 
lization; 
evacuation  is 
neglected 

Full 
therma¬ 
lizat  ion; 
basic 
evacu¬ 
ation 

Full  thermal! za- 
tion;  full  evacu¬ 
ation 

A 

2 

1  | 

1.412 

1 _ i 

j  1.047 

0.835 

0.696 

A  practical  example  for  the  double-pass  amplification  is 
given  in  the  third  section  of  this  article. 

We  will  study  the  amplification  problem  for  pulse  series 
in  the  following.  The  physical  process  and  calculation  method 
are  basically  no  different  from  the  double-pass  amplification. 
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However,  the  number  of  pulses  is  much  greater;  hence,  the 
expression  is  more  tedious  for  calculating  the  gain  of  every 
pulse.  Accordingly,  the  optimum  expression  will  be  the  recur¬ 
rence  relation.  Besides,  the  interval  between  pulses  is  usually 
shorter  ;  hence,  the  evacuation  effect  of  the 

final  state  may  not  be  considered.  Consequent ly ,  the  variation 
of  particle  numbers  must  be  described  by  three  quantities.  The 
recurrence  relation  Equation  (8)  is  reduced  to: 

T+  _ £ul!L 

L  ,L  -  j  X'2  i  *r  ki 

X(1 + 

The  following  initial  conditions  are  given: 


~  e' »?» ■=  0,  s'P  -  (1  +  ka)  Np, 

(k) 

The  small  signal  gain,  g  ,  and  energy  density  of  output 
(  r  )  ° 

light,  E  (L),  of  any  pulse  can  be  derived,  leading  to: 

E(t> (L)  -  E. In  [1  +  (eF?VB*  - 1)  <tf°]  e~“L0 

The  above  formulas  apply  to  any  case  which  satisfies 
5/r«cl  .  It  is  not  difficult  to  prove  that  the  above  for¬ 
mulas  are  reduced  to  the  results  derived  in  Ref.  2  for  the  special 
case  of  5/T»l  and  O„o 

II.  THE  ENERGY  EXTRACTION  EFFICIENCY  OF  A  HIGH  POWER  AMPLIFIER 

Energy  conversion  efficiency  Is  an  important  characteristic 
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indicator  of  an  amplifier.  We  will  not  discuss  here  the  problem 
related  to  light  pumping  efficiency,  but  we  will  study  the 
efficiency  of  storage  energy  of  initial  reversed  particles  converted 
to  laser  energy  —  extracting  energy.  The  so-called  "extraction 
laser  energy"  means  the  difference  between  the  output  and  input 
light  energies  of  an  amplifier.  It  equals  the  light  energy- 
released  through  excited  transitions  after  subtracting  the  energy 
loss  of  operating  materials.  For  convenience,  we  define  a  local 
extraction  efficiency  n(z)  and  a  total  extraction  efficiency  n: 

/  _  extraction  energy  density  at  point  Z  in  the  ajrrdifler 

storage  energy  density  at  point  Z 


hv  (fl.-w  —  —  aE 

.  hvN0 

total  extraction  energy  in  an  amplifier 
total  storage  energy 


(13) 


-  \\  (El  -  Eo)  dS ! j|  ^hvNodZ)  dS  {lh) 

•  « 

where  n^Q  is  the  density  of  the  initial  particles  in  the  metasta¬ 
ble  state;  n^  is  the  density  of  remaining  particles  in  the  metas¬ 
table  state  after  the  laser  beam  passes;  Nq  is  the  density  of 
initial  reversed  particles;  Er  ,  E  are  densities  of  output  and 
input  light  energies,  respectively. 

The  total  extracting  efficiency  can  also  be  rewritten  in 
the  form; 


v~\\(]]-n{E)N0{Z)dz)dS 

a 

m:  N0(Z)dZ  'jdS 

a 


It  can  be  seen  that  n  is  the  mean  value  of  »)(£) 


We  first  discuss  the  local  extraction 
n(Z)  can  be  expressed,  using  Equation  (1),  as 


efficiency  problem, 
the  following: 


2 


1  a  E(Z ) 

2  /30 


(15) 


where  F'.(Z)  is  a  variable,  is  a  parameter,  and  the  variation 

rule  of  n(Z)  is  shown  in  Figure  3-  The  curves  in  the  figure  reflect 

the  fundamental  characteristics  of  extraction  efficiency.  As  can 

be  seen,  for  any  given  energy  density,  extraction  efficiency  increase 

with  increasing  /30/<*  .  As  a  result,  the  intense  light  pumping 

is  in  general  advantageous  for  upgrading  efficiency.  Secondly, 

for  a  given  &/<*.  ,  there  exists  a  maximum  on  each  curve.  Thus,  in 

order  to  achieve  a  high  efficiency,  the  operating  region  in  an 

amplifier  (i.  e.  the  variation  range  from  input  energy  density,  E  , 

to  output  energy  density,  EL)  must  be  located  near  the  maximum. 

The  efficiency  can  be  reduced  when  light  energy  density  is  too 

large  or  too  small.  The  light  energy  density  corresponding  to  the 

maximum,  E  ,  can  easily  be  derived: 
irm 


^.-jJ.lnC/So/^o 


The  value  of  A> /«  mostly  observed  in  experiments  is  around  10. 

The  corresponding  E  is  over  double  the  saturated  energy  density 
E  ,  which  far  exceeds  the  loading  capability  of  neodymium  glass. 
Thus,  the  condition  for  optimum  operating  region  is  hardly 
realized.  In  the  following,  we  will  study  a  problem  which  has 
a  practical  significance:  how  to  increase  extraction  efficiency 
without  increasing  light  energy  density.  Assuming  that  the 
transverse  distribution  is  uniform,  the  total  extraction  efficiency 
can  be  reduced  to: 


(l?t  —  Eo) /2E,/3()L 


(16) 


J)  JE(Z)(R4/H*,> 


Figure  3.  Local  extraction  efficiency. 
1)  E(Z)  (joules/cm^) 


There  are  four  variables  affecting  n:  E^,  Eq,  B^,  L,  but  only 
three  of  them  are  independent  through  Equation  (1).  To  be  ex¬ 
plicit,  Equation  (1)  can  be  transformed  into  the  following  expres¬ 
sion  : 


L  ~  (*‘  dE/  [i?,/90  (1  -  •-»*•)  -  aft) , 

}*• 

which  is  substituted  into  Equation  (16)  to  result  in: 


(17) 


E0)  /2  E,  (/So/a)  dE 

/[E,  (fio/ a)  (1  —  E  *•)  —  Ty )  0 


(18) 


When  the  operating  region  is  determined,  n  varies  with  f}0/a 
It  can  be  proved  that  only  if  Et,>E0 — fih--->0o  is  it  always 

_  ’  d  (po/ a) 

effective.  That  is,  n  increases  with  increasing  p0/a  .  When 
Po/a~¥<x.  >  v  approaches  a  maximum  7im  : 


Vm  -  (El  -  E0)  /2E,1n  [  (•**'■•-1) 


(19) 


Several  sets  of  relation  curves  are  shown  in  Figure  M , 

calculated  for  several  popular  experimental  environments.  Those 
curves  indicate  that  an  upgrading  of  /80/a  will  enhance  the  effi¬ 
ciency.  It  is  particularly  prominent  when  fo/a  is  small.  But 
when  /So/«  approaches  10,  it  will  be  close  to  the  maximum  extraction 
efficiency,  and  any  further  increase  in  A>/«  cannot  help  the 
efficiency  significantly.  Consequently,  how  to  choose  &,/<* 
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value  properly  Is  a  critical  problem  needed  to  be  taken  into 
account  in  amplifier  design. 

When  the  /9 „/«*  is  given,  the  length  of  the  amplifier  cannot 

be  random  in  order  to  satisfy  the  input  and  output  requirements. 

The  length  is  determined  by  Equation  (27).  It  can  easily  be  seen 
that  L  decreases  with  increasing  S  .  Subject  to  the  limitation  of 
amplifier  loading,  we  should  adopt  intense  light  pumping  (large 
and  short  bar  (small  L))  along  with  a  value  of  & Ja  no 
greater  than  10,  in  order  to  achieve  a  high  efficiency. 

III.  TWO  ROUND,  THREE  PASS  DISK  AMPLIFIER 

We  designed  and  developed  a  large-scale  disk  amplifier  as 
an  application  of  the  foregoing  theory.  The  amplifier  is  composed 
of  six  slices  of  neodymium  glass;  each  slice  has  a  thickness  of 
35  mm,  width  of  180  mm,  height  of  300  mm.  The  slice  is  inclined 
at  a  Bragg's  angle  (56°24').  The  dimensions  of  the  light  pass 
aperture  are  100  mm  x  300  mm.  The  amplifier  also  includes  12 
upright  xenon  lamp  light  pumps  as  shown  in  Figure  5.  The  light 
beam  passes  the  lower,  middle  and  upper  parts  consecutively,  thus 


Figure  4.  The  relation  between  total  extraction 
efficiency  and  <V«»  .  Key:  1)  operating  region. 
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forming  three  passes.  In  order  to  raise  the  gain  and  efficiency, 
another  passage  is  repeated.  Figure  5(B)  shows  the  light  paths 
following  the  principle  of  two  round,  three  pass  operation.  In 
order  to  fully  utilize  the  advantages  of  double-pass  amplification 
and  final  state  evacuation  effect,  the  time  interval  between  the 
two  consecutive  rounds  of  amplification  must  be  as  large  as  pos¬ 
sible.  Consequently,  the  distance  between  reflection  mirrors  is 
about  15  m,  and  $»t,  is  taken  for  calculation  of  the  gain  in 
the  second  round  of  amplification. 


Figure  5.  Disk  amplifier:  A)  structure; 

B)  conceptual  light  paths;  1)  xenon  lamp; 

2)  light  concentrator;  )  neodymium  glass 
slice;  4)  disk  amplifier;  5)  input;  6)  output; 

7)  15  m. 

The  static  absorption  coefficient  of  the  neodymium  glass 
is  ao*»0.16%  .  However,  for  light  pumping  operation,  the  dynamic 

absorption  for  the  1.06  ym  laser  is  larger  than  °u  (about  two 
to  three  times  a„  )  because  of  the  upward  excited  transitions 
of  instantaneous  color  centers  and  metastable  states,  as  well  as  the 
nonlinear  absorption  during  intense  laser  interaction.  We  take  the 
absorption  coefficient  to  be  “  2. ^-0.4%, 

Surface  loss  is  determined  by  experiments.  The  total  loss  of  each 
neodymium  glass  slice  thus  obtained  is  2%. 

The  Initial  gain  coefficient  is  derived  from  small-scale 

simulated  experimental  data.  When  the  density  of  light  pumping 

3 

(electric  energy /neodymium  glass  volume)  is  29  joules/cm  ,  6Q  = 

2 . QG%  cm  ^ . 

14 


In  accordance  with  formulas  previously  derived,  we 
calculate,  for  a  typical  case  with  input  laser  energy  of  40 
joules,  the  gain  and  output  energy  of  each  round  of  amplifica¬ 
tion  which  is  listed  in  the  following  table: 


Energy 

Gain 

Output 

Energy 

Average  Residual 
at 

Reversed  Number 

Footnotes 

First  round 

4.24 

times 

172  joules 

CO 

U1 

Diameters  of  light 
beam  are  4>70,  80, 
and  90,  respectively 

Second  round 

2.50 

430 

64$ 

Diameter  of  light 
beam  is  4>100 

Total  of 
two  rounds 

10.6 

430 

64$ 

Reversed 

later 

3.01 

127 

Reflectivity  of 
target  surface 
is  10$ 

The  experimental  results  are:  input  energy  40  joules, 
output  energy  399  joules,  gain  10.  Experimental  values  and  values 
obtained  by  theoretical  calculations  are  basically  in  agreement. 

The  disk  amplifier  is  intended  for  use  in  a  nuclear  fusion 
experiment,  hence  the  table  also  lists  some  calculated  results  for 
the  gain  of  the  reversed  laser.  In  the  calculation,  we  assume  the 
retardation  rate  of  each  pass  as  8$.  Experimental  values  are  in 
agreement  with  calculated  values. 
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In  order  to  further  examine  the  theory,  we  have  conducted 
measurements  and  theoretical  calculations  of  gains  for  various 
energies  of  light  pumping  and  various  laser  intensities.  The 
results  are  listed  in  the  following  table.* 


Series  ft 

Density  of 
Light  Puiping 

Input 

Laser  Energy 

Gain 

Coefficient 

Calculated 
Values 
for  Gain 

Experi¬ 

mental 

Values 

1 

O 

11.4  joules/cm 

44.5  joules 

1.37%  cm-1 

2.30 

times 

2.80 

times 

2 

11.4 

50.5 

1.37* 

2.23 

2.76 

3 

16.4 

62.0 

1.88* 

3.64 

3-98 

4 

23.6 

53-5 

2.54^ 

7.05 

7.00 

5 

26.2 

44.0 

2.71 ! 

9.06 

8.40 

6 

29.0 

40.0 

2.86 % 

10.6 

9.97 

The  results  indicate  that  under  various  operating  conditions 
the  calculated  values  are  basically  in  agreement  with  experiments. 

Finally,  we  want  to  examine  the  problem  with  respect  to  the 
effic  ency  of  energy  conversion.  When  the  density  of  light  pumping 
is  29  joules/cm'),  the  total  input  electric  energy  is  3-6  x  lCr 
Joules,  and  the  total  extraction  laser  energy  in  two  rounds  of  ampli¬ 
fications  is  360  joules,  thus  leading  to  a  total  efficiency  of  0.1". 
The  efficiency  of  electric  energy  converted  to  the  storage  energy  of 

■* - 

NOTE:  The  pain  coefficient  is  not  simply  proportional  to  the  light  punping 
energy.  The  effect  due  to  the  variation  of  light  emission  of  a  xenon  lanp  must 
be  taken  into  account . 


reversed  particles  (i.e.  the  efficiency  of  light  pumping)  is  0.89%, 
and  the  extraction  efficiency  is  11.2%.  Compared  to  other  similar 
types  of  amplifiers,  the  extraction  efficiency  is  higher  for  this 
case  due  to  the  use  of  the  double-pass  amplification  along  with  a 
higher  energy  density  of  the  light  beam.  The  C  series  disk  installed 
in  the  Cyclops  at  Livermore  Laboratory  in  the  United  States  has 
a  light  pumping  efficiency  of  close  to  1%.  However,  the  energy 
density  of  the  light  beam  can  be  controlled  only  in  the  range  of 
a  small  signal,  due  to  the  restriction  of  the  nonlinear  effect. 

Thus,  the  extracting  efficiency  is  very  low,  only  about  1% .  It 
is  expected  that,  if  our  device  is  applied  to  the  amplification 
of  pic  osecond  pulses,  the  efficiency  will  decline. 

The  foregoing  discussions  imply  that,  in  order  to  achieve 
higher  efficiency,  research  on  high-power  amplifiers  must  result  in 
the  improvement  of  light  beam  transmission  conditions  and  over¬ 
coming  the  nonlinear  restriction  so  as  to  raise  the  loading  capa¬ 
bility  of  operating  materials  as  much  as  possible.  It  is  wise  to 
emphasize  light  pumping. 
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Abstract 


Based  on  practical  work  in  1971-1976,  some  problems  to  be 
considered  in  designing  of  rapid  expansion  nozzles  for  gas-dynamic 
lasers  were  summarized,  including  the  selection  of  important  param¬ 
eters,  the  discussion  on  methods  for  contoured  nozzles  design  and 
the  correction  for  boundary  layer  growth,  etc.  Through  a  large 
amount  of  numerical  analysis,  a  simple  convenient  method  for  cor¬ 
recting  boundary  layer  growth  under  various  conditions  was  also  pre¬ 
sented  . 

I.  INTRODUCTION 

The  rapid  expansion  nozzle  is  a  key  component  for  the 
gas-dynamic  laser.  Its  major  function  is  to  rapidly  and  effec¬ 
tively  freeze  the  vibration  energy  in  gas-dynamic  operation.  It 
is  different  from  the  conventional  nozzle  installed  in  a  super¬ 
sonic  wind  tunnel.  The  basic  principle  in  designing  a  gas-dynamic 
nozzle  in  a  laser  device  should  be:  only  the  premise  guaranteeing 
an  air  current  with  minimum  perturbation.  Maximum  freezing  should 
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be  achieved.  The  perturbation  of  an  air  current,  i.e.  the  non¬ 
uniformity  of  an  air  current,  will  directly  affect  the  quality  of  the 
output  light  beam,  while  the  freezing  effect  is  related  to  the 
efficiency  of  the  device.  There  is  some  contradiction  between 
guaranteeing  a  uniform  air  current  and  effective  freezing.  It 
requires  overall  balance  and  reasonable  adjustment. 


II.  FACTORS  WHICH  AFFECT  THE  FREEZING  EFFECTS  OF  A  NOZZLE 

A  large  amount  of  theoretical  analysis  and  experimental 
studies  have  indicated  [1-8]  that,  for  a  certain  air  current  re¬ 
tardation  condition,  the  major  factor  which  affects  the  freezing 
is  the  geometry  of  the  nozzle.  Its  primary  geometrical  parameters 
include  curvature  radii  of  the  throat,  R^  and  ,  the  maximum 
expansion  angle  of  the  air  current,  0^,  the  height  of  the  throat 
pipe,  R*,  and  the  area  ratio,  A^/A*  (see  Figure  1). 

How  the  shape  of  a  subsonic  nozzle  affects  the  freezing 
effect  cannot  be  overlooked.  A  small  curvature  radius,  R^ ,  of 
the  throat  and  a  large  entrance  angle,  0,,  for  an  air  current  enable 
the  air  current  to  become  unbalanced  in  the  subsonic  range  and  this 
expedites  the  freezing  of  vibration  energy.  R^=0  and  0^=90°  repre¬ 
sent  the  fastest  subsonic  expansion  and  the  best  freezing  effect. 

But  at  R-L=h*  and  0^=45°,  the  degradation  of  quality  is  not  very 
significant  [6].  The  curvature  radius,  R^,  of  the  throat  in  the 
supersonic  portion  of  a  nozzle  has  the  same  influence  on  the  freezin 
effect  as  R^ .  In  principle,  the  smaller  the  better.  Its  value  is 
determined  by  the  calculation  method  for  the  contoured  nozzle.  Ir. 
general,  the  minimum  length  of  a  nozzle  is  calculated  by  the  method 
of  characteristic  lines  with  which  the  throat  portion  appears  at  an 
acute  angle  and  R2=0.  If  the  Foelsch  method  is  applied  to  calcu¬ 
lating  the  nozzle  [9,10],  the  curvature  radius  of  the  throat  will 
be  related  to  h*  and  G^ .  In  general,  lia<h *  .  As  for  3^,  the 

maximum  value  of  0^  allowed  in  the  theory  is  generally  taken,  i.e. 
half  of  Prandt  1-Moyer  expansion  angle  01  £  =  in  designing  Mach 

correspondence . 
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Figure  1.  Geometrical  parameters  of  a  nozzle. 

The  height  h*  of  the  nozzle  throat  is  an  important  parameter 
which  characterizes  the  freezing  efficiency  of  the  nozzle.  For 
nozzles  with  the  same  area  ratio  (viscosity  is  not  taken  into  account), 
the  geometries  are  similar.  The  smaller  the  dimensions  of  the 
throat  pipe,  the  shorter  the  length  of  the  nozzle,  and  the  varia¬ 
tions  of  air  current  parameters  along  the  axial  direction  are 
larger.  Consequently,  the  freezing  effect  is  more  prominent.  Small 
throat  height  is  a  characteristic  of  a  nozzle  of  a  gas-dynamic  laser 
device,  usually  on  the  order  of  0.5  mm.  If  manufacture  and  assembly 
can  be  warranted,  and  if  the  contour  surface  of  the  nozzle  can  be 
sufficiently  cooled  (especially  in  the  throat  region),  the  height 
of  the  throat  can  be  further  reduced,  e.  g.  reduced  to  as  small  as 
0.1  mm  [11]. 

The  ratio  of  the  nozzle  outlet  area  to  the  throat  area  AM*, 
is  another  important  parameter  for  determining  the  freezing  efficiency 
of  a  nozzle.  The  larger  the  area  ratio,  the  larger  the  degree  of 
air  current  expansion,  and  the  lower  the  temperature  and  pressure 
of  the  air  current  in  the  nozzle.  In  general,  the  area  ratio  is 
determined  by  the  relation  of  equal  entropy  in  accordance  with  the 
ratio  of  retarding  temperature,  T  ,  accessible  in  the  burning 
chamber  to  regulated  temperature  of  the  air  current  in  a  light 
cavity,  T.  The  air  current  temperature  in  a  light  cavity  is 
regulated  to  guarantee  the  vapor  will  not  be  condensed.  In  the 
past ,  T:s300oK  was  usually  taken.  But  the  vapor  in  a  nozzle  is  not 
condensed  immediately  after  saturation  during  the  process  of  rapid 
expansion.  It  generally  requires  i40"80oC  excessive  cooling  tem¬ 
perature.  In  addition,  part  of  vibration  energy  in  CO^-H^  is 
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finally  converted  to  kinetic  energy  through  a  relaxation  process, 
leading  to  a  higher  air  current  temperature  compared  to  the  calculated 
temperature  with  the  equal-  entropy  relation.  Accordingly,  by 
considering  these  two  factors,  the  area  ratio  can  be  chosen  larger. 

In  the  discussion  of  factors  which  affect  the  freezing 
effect,  the  problem  regarding  the  recombination  parameter  pQh* 
must  be  taken  into  account.  Stollery  and  Smith  [12]  proposed  a 
similar  parameter  PQh*  to  describe  the  freezing  effect  before, 
when  they  studied  the  non-equilibrium  effect  of  a  nozzle.  Gerry [1] 
derived  the  fact  that  freezing  usually  occurs  near  the  throat  region. 
He  derived  a  principle  which  describes  the  freezing  effect  of  a 
gas-dynamic  laser  device,  represented  by  pQh*: 

jj0A*<2(pr<1#)!P.u* 

This  principle  was  supported  by  Anderson  [ 5] .  He  deliberated  the 
implication  of  this  recombination  parameter  and  pointed  out  that  if 
other  parameters  of  the  nozzle  (such  as  gas  components,  temperature, 
area  ratio,  etc.)  are  invariant,  the  maximum  gain  and  the  maximum 
usable  energy  are  all  single-valued  functions  of  pQh*.  He  also 
suggested  that  it  is  more  appropriate  to  take  pQh*~l  atm.  cm.  We 
believe  that  this  recombination  parameter  can  be  used  as  a  reference 
In  designing. 

m .  FACTOI.S  WHICH  AFFECT  THE  UNIFORMITY  OF  THE  AIR  CURRENT 

There  are  many  factors  which  affect  the  uniformity  of 
the  air  current.  For  arrayed  nozzles  of  a  gas-dynamic  laser 
device,  major  factors  are  wave  series  and  tail  trace.  The  tail 
of  thin  disks  in  the  arrayed  nozzles  must  be  cut  due  to  considera¬ 
tion  of  structure  strength,  thus  inevitably  generating  excited 
waves  of  certain  intensity.  The  air  current  through  the  tail  will 


generate  fluctuating  tall  traces  due  to  viscosity.  Secondly,  im¬ 
proper  design  of  a  nozzle  (including  correction  of  the  viscosity 
effect)  and  improper  assembly  may  induce  air  currents  to  generate 
various  perturbations  (air  current  separation,  various  degrees  of 
expansion  and  compressed  wave  series,  as  well  as  the  interaction 
between  the  wave  series  and  boundary  layer).  In  addition,  the  non¬ 
uniformity  and  instability  of  the  initial  current  from  the  burning 
chamber,  due  to  burning  and  mixing,  must  be  taken  into  account. 

Quantitative  evaluation  of  the  effect  of  these  factors  on  the 
uniformity  of  air  current  is  a  difficult  and  complicated  problem. 

To  evaluate  the  effect  of  a  non-uniformity  of  the  air  current 
incurred  by  these  factors  is  even  more  difficult .  We  can  only 
offer  a  crude  qualitative  introduction. 

Simons[l8]  pointed  out  that,  if  the  boundary  layer  on  the 
side  wall  of  the  nozzle  is  not  corrected,  the  excited  waves  induced 
by  this  type  of  boundary  layer  may  make  the  air  current  generate  a 
10%  density  non-uniformity.  A  simple  approach  which  enables  side 
walls  to  expand  outward  linearly  (precise  quantity  is  determined 
by  calculation,  about  on  the  order  of  0.5  ~  1°),  can  reduce  such 
perturbations  by  more  than  half.  Clark[l4]  estimated,  with  the 
method  of  geometrical  optics,  the  effect  of  such  perturbation  on 
the  quality  of  a  light  beam  for  a  co-focus  unstable  cavity.  When 
Jp/p- 7%  ,  remote  fields  parallel  to  and  perpendicular  to  the  air 

current  are  greatly  deformed.  In  an  analysis  of  this  problem, 
Russell[15]  pointed  out  that  the  density  perturbation  induced  by 
excited  waves  near  side  walls  must  be  limited  to  less  than  1%. 

For  an  excited  wave  intersecting  the  optical  axis,  it  can 
only  deflect  light  beams  in  general,  i.e.  the  phase  deviation  is 
linear.  The  intensity  of  its  remote  field  is  unchanged  relative 
to  the  deflected  optical  axis.  Accordingly,  such  a  deviation  can 
be  corrected  through  an  external  light  path  system.  However,  these 
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excited  waves  interact  with  tail  trace  or  boundary  layer,  pos¬ 
sibly  leading  to  a  high-order  phase  deviation.  Therefore,  the 
intensity  of  this  type  of  excited  waves  must  be  reduced  as  much 
as  possible. 


Since  the  optical  axis  of  gas-dynamic  laser  devices  is 
arranged  perpendicular  to  the  rear  part  of  the  nozzle  disk,  the 
tail  trace  and  fluctuating  current  do  not  affect  the  quality  of 
the  light  beam  very  significantly [ 14 ] .  They  affect  the  quality  of 
the  light  beam  primarily  through  the  divergency  of  the  light  beam. 

To  Properly  evaluate  these  effects  is  very  complicated. 

The  manufacture  of  nozzles,  particularly  those  of  the 
throat  portion,  must  be  accurate  and  provide  excellent  surface 
smoothness.  Wagner[l6]  suggested  that  the  accuracy  in  manufac¬ 
turing  must  be  guaranteed  to  -  0.025  mm,  and  the  surface  smooth- 
ness  must  be  0.3  x  10~  m  (corresponding  to  711).  It  can  only  be 
achieved  through  delicate  grinding  and  projection.  Currently  we 
do  not  have  technology  to  achieve  that  goal. 

The  assembly  of  nozzles,  especially  accurate  dimensions  of 
the  throat  portion,  is  very  important.  Sin^e  the  dimensions  of 

-f. 

the  throat  portion  are  very  small,  a  small  assembly  error  (-  0.02) 
may  lead  to  large  relative  error  (-  10$),  hence  seriously  affecting 
the  uniformity  of  the  air  current. 

In  order  to  have  a  uniform  current,  we  must  also  assure  the 

stability  of  nozzle  dimensions  during  operation.  As  a  result, 

nozzle  cooling  is  a  necessary  step.  In  addition,  because  of  the 

viscosity  effect,  the  air  current  near  the  surface  of  the  nozzle 

wall  will  be  decelerated,  converting  the  kinetic  energy  into  thermal 

energy.  The  air  current  will  be  thereby  heated  and  there  may  be 

an  unfavorable  impact  on  the  reversal  of  particles.  For  a  nozzle 

with  To  =  1^00  K,  PQ =2* 5  atm,  M=4 ,  h*=0.2  mm,  the  loss  of  gain  may 

reach  21$[ll]  when  the  temperature  at  the  wall  surface  Is  T  =700K. 

w 
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Thus,  from  the  viewpoint  of  freezing,  the  nozzle  must  be  suf¬ 
ficiently  cooled. 


IV.  DISCUSSIONS  ON  THE  DESIGN  OF  NOZZLES 

In  order  to  reduce  perturbations,  a  detailed  gas-dynamical 
calculation  of  the  nozzle  wall  is  necessary.  The  guidelines  in 
treatment  of  this  problem  may  follow  that  of  a  nozzle  installed  in 
a  conventional  supersonic  wind  tunnel.  Since  the  vibration  energy 
is  only  a  few  percents  of  the  total  enthalpy,  when  the  ratio  of 
gas  specific  heats,  y  ,  is  given,  the  non-equilibrium  effect 
can  be  neglected.  The  gas-dynamic  parameters  are  not  affected 
significantly  if  nozzle  calculations  are  performed  based  on  ideal 
air  currents  of  equal  entropy[ll]. 

The  design  for  the  subsonic  part  is  much  easier.  Only  a 
linear  contraction  of  nozzle  areas  and  no  counter-pressure  in 
the  axial  direction  need  be  assured.  It  is  not  necessary  to 
calculate  the  wall  shape  in  detail[15].  For  example,  even  for  the 
extreme  case  with  entrance  angle  of  0^=90°  and  a  curvature  radius  of 
the  throat,  R^=  0,  air  streams  separation  is  not  observed,  and  the 
requirement  for  the  uniformity  of  the  air  stream  is  sat isf ied[ G ] . 

In  addition,  the  air  current  subsonic  part  is  very  close  to  a 
one-dimensional  current,  only  if  0.  Therefore,  the  linear  con¬ 

traction  or  wedge  contraction,  i.e.  6^=45°,  R-^-h*,  are  normally 
applied  to  the  design  of  the  subsonic  unit. 

For  nozzles  in  the  supersonic  portion,  reportedly  the  nozzle 
with  minimum  length  calculated  by  the  characteristic  line  method 
has  been  adopted  in  all  cases  [5>  6,  15,  16].  In  the  design,  the 
throat  portion  is  assumed  to  be  a  straight  sonic  line,  and  the 
initial  expansion  angle  is  taken  as  the  maximum  value,  o ^  = 

The  throat  pipe  of  this  type  of  nozzle  appears  at  an  acute  angle, 
along  with  maximum  expansion  speed  and  minimum  length.  In  terms 
of  freezing  vibration  energy,  it  is  certainly  a  better  design. 
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But  there  are  problems  in  practice.  Firstly,  the  requirements  in 
manufacture  and  assembly  are  very  strict.  For  an  acute  angle  throat 
pipe,  if  the  dislocation  of  the  axial  direction  3s  lj  of  throat 
height,  a  corresponding  rotation  angle  perturbation  of  1°  in  the 
air  current  is  generatedf 15 ] .  Secondly,  the  acute  angle  may 
easily  cause  separation  of  the  air  current,  leading  to  rapid  growth 
of  the  boundary  layer  on  the  side  wall  in  the  throat  portion  and 
generation  of  excited  waves  in  the  threat.  The  assumption  of  a 
straight  sonic  line  also  leads  to  over-expansion  in  the  lower 
portion  of  the  throat  pipe,  which  may  cause  separation  and  generate 
excited  waves.  Accordingly,  Anderson  suggested  using  the  arc 
shape  to  replace  the  acute  angle.  This  approach  may  reduce  theoreti¬ 
cal  gain  and  maximum  usable  energy  by  approximately  15%,  but  the 
real  benefit  may  surpass  it. 

Consequently,  the  nozzle  of  minimum  length  designed  by  the 
characteristic  line  method  is  not  a  best  one.  Vie  suggest  that  a 
method  (Foelsch  method)  be  employed  which  has  been  most  popular  in 
the  design  of  conventional  nozzles  installed  in  a  supersonic  wind 
tunnel  in  China.  The  method  assumes  that  the  air  current  before 
the  nozzle  turning  point  has  a  spring  current-like  distribution. 

It  is  a  very  convenient  analytical  method  which  applies  the  principle 
of  characteristic  lines  to  calculate  the  nozzle  well  [10]  as 
shown  in  Figure  2.  Major  formulas  for  calculation  are  listed  below. 


Figure  2.  Design  diagram  of  nozzle. 
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The  nozzle  designed  by  the  Foelsch  method  has  a  curvature 
radius  in  the  throat  portion  of  less  than  the  throat  height.  The 
method  takes  care  of  requirements  for  freezing  and  uniformity  of 
the  air  current.  Besides,  its  wall  coordinates  are  represented  by 
simple  analytical  expressions  which  are  convenient  for  various 
computations . and  analyses.  This  method  provides  adaptability  of 
gas-dynamic  laser  nozzles.  Including  higher  Mach  number  and  a 
larger  expansion  angle.  Comparisons  have  been  made  with  the 
characteristic  line  method  by  using  electronic  computers.  Results 
indicate  that,  even  though  the  spring  current  assumption  has  some 
discrepancies  under  this  condition,  the  method  is  still  usable. 
Results  for  practical  operation  of  several  types  of  gas-dynamic 
laser  devices  have  indicated  the  following:  the  small  signal  gain 
measured  has  reached  or  even  surpassed  the  best  values  of  similar 
devices  in  other  countries. 
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For  gas-dynamic  laser  devices  operated  at  low  pressure,  the 
density  perturbation  is  not  relatively  serious  because  of  the  low 
cavity  pressure.  Thus  the  nozzle  wall  need  not  be  carefully  con¬ 
sidered.  Nozzles  in  simple  wedge  shape  or  cone  shape  can  be 
used  [15]. 

V.  CALCULATIONS  FOR  FREEZING  > 

For  a  gas-dynamic  calculation,  y  is  a  basic  parameter. 

There  are  two  viewpoints  on  Y  •  One  regards  that  y  does  not 
affect  the  theoretical  wall  significantly  and  Y  =1.4  can  be 
assumed  in  designs.  Its  basic  rationale  is:  once  area  ratios 
are  identical,  walls  calculated  by  the  same  method  but  with 
different  y  will  not  be  much  different .  Another  viewpoint 
is  completely  different.  It  considers  y  value  critical  in 
wall  calculations  [15]. 

We  believe  that  taking  y  as  a  constant  is  an  approxi¬ 
mation,  but  we  realize  that  Y  affects  the  Mach  number  of  the 
nozzle  and  the  displacement  thickness  of  boundary  layer  in  some 
degree.  On  the  other  hand,  it  is  inappropriate  to  choose  a  Y 
value  randomly.  Besides,  the  fact  that  nozzle  wall  coordinates  are 
not  significantly  affected  by  the  variation  of  y  may  provide  a 
basis  for  experiments  to  search  for  optimum  parameters  when  gas 
components  are  changed  in  a  fixed  nozzle.  However,  such  cases 
with  optimum  parameters  usually  do  not  apply  to  the  case  of  nozzle 
design.  In  practical  calculations,  conceptual  "freezing  y  "  may 
apply,  l.e.  assume  that  mixed  air  current  is  an  equilibrium  current 
before  it  flows  through  the  throat  pipe,  while  gas  vibration  energy 
is  suddenly  frozen  in  the  throat  pipe. 

In  consideration  of  the  high  characteristic  temperature  of 
two  vibration  energy  levels  for  H2O,  the  vibration  energy 

for  these  two  enerj^/  levels  is  negligible  for  a  given  retardation  condition. 


Thus,  starting  from  the  throat  portion,  only  vibration  energies 
of  energy  levels  in  CO^  and  level  in  I^O  vary  with  the 

temperature  of  the  air  current  in  a  mixed  gas. 
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Since  E  (T)  is  a  function  of  temperature,  the  specific  heat 
at  constant  volume  Cv  of  the  mixed  gas  (and  thus  y  )  becomes  a 
function  of  temperature.  Accordingly,  if  y  is  taken  constant, 
we  must  find  the  average  value  of  in  some  sense. 

Calculation  steps  are  as  follows: 


1)  Following  the  approach  of  equi-entropy  variable  specific 
heat,  calculate  the  temperature  of  the  air  current,  T# , 
in  the  nozzle  throat,  which  is  used  as  the  freezing 
temperature  of  upper  energy  levels  of  mixed  gas,  and  is 
the  initial  temperature  used  to  calculate  the  freezing  y 


2)  Calculate  the  average  value  of  Cv  by  applying  the 


mean  value  theorem, 


(Cv),“2V^Ir!  {°')tdT 

Ev(T.)-Et(T.) 
T.-T, 


(*5) 


where  Tg  is  the  air  current  temperature  at  the  outlet 
of  the  nozzle. 

3)  Freezing  7  can  be  calculated  by  the  relation, 

Cv  "  (Cv}Tr  +  (5vV 


7-1  + 


rt_ 

Q. 


(6) 


28 


where  (COjv—tt  -HO/mo.  +  '/'(O  +  &R'l‘n&,  ,  R  is 

0 

the  universal  gas  constant,  «^r(  is  the  component. 

VI.  CALCULATIONS  OF  NOZZLE  BOUNDARY  LAYER 

For  rapid  expansion  nozzles  in  gas-dynamic  laser  devices, 
the  development  of  the  boundary  layer  in  the  axial  direction  of 
the  nozzle  is  very  fast  because  of  the  very  small  dimensions.  As 
a  result,  the  air  current  uniformity  is  substantially  affected  by 
viscosity . 

Whether  the  boundary  layer  in  the  nozzle  of  a  gas-dynamic 
laser  device  is  laminar  or  turbulent  is  a  problem  bearing  further 
discussion.  In  terms  of  the  Reynolds  number  of  a  typical  gas- 
dynamic  laser  nozzle  (nozzle  length  is  the  characteristic  length 
in  this  case),  which  is  generally  smaller  than  10^,  it  should  be 
laminar  in  principle.  But,  in  fact,  because  of  very  small  curvature 
radius  of  the  throat  part  (acute  angle)  which  leads  to  large  expan¬ 
sion  angle  for  air  current,  the  air  current  may  separate  and  become 
a  turbulent  boundary  layer.  Foreign  publications  in  this  regard  are 
mostly  based  on  the  laminar  calculations  [12,13,15,17].  However, 
early  experiments  conducted  by  the  NOL  in  the  United  States  indicated 
that  the  boundary  layer  in  the  nozzle  is  turbulent. 

The  Cohen  and  Reshotko[17]  approximation,  which  has  been 

very  extensively  used  in  boundary  layer  calculations  in  gas-dynamic 

laser  devices  in  foreign  countries,  will  be  adapted.  We  separate 

it  into  two  cases:  heat  conduction  wall  (t  /t  =0.4)  and  insulated 

w  o 

wall  (t  /t  =  1).  Three  different  ratios  of  specific  heat  (  1  =1.33, 

w  o 

1.36,  1.4)  and  five  different  Mach  numbers  (M-4.0,  4.25.  4.5,  4.75,  6.0).  are 
used  to  calculate  parameters  sucn  as  displacement  tnickness 5*  for  a  boundary  layer 
in  a  nozzle  which  has  been  calculated  by  the  Foelsch  method. 


Even  nozzles  with  similar  geometrical  structure  may  develop 
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different  boundary  layers.  Therefore,  the  displacement  thick¬ 
ness  for  a  nozzle  boundary  layer  cannot  be  considered  directly 
as  a  simplified  similar  parameter,  but  only  can  be  expressed  in 
a  combination  parameter,  such  as  CR»)  '  » 

For  convenience,  we  adopt  a  combination  parameter  for  the  dis¬ 
placement  thickness  in  the  form: 


DELTAFI1T 


7,  Mb, 


(7) 


Once  the  combination  parameter  for  displacement  thickness 
is  given,  we  can  find  the  displacement  thickness  along  the  nozzle 
wall  surface,  8*(s)  ,  based  on  various  given  parameters,  such  as 

7,  -Me,  h\  T0>  po,  <«/A>  .  See  Figure  3- 


Figure  3.  Boundary  layer  correction:  1)  wall 
after  correction;  2)  theoretical  wall. 


The  coordinates  of  the  nozzle  wall  after  correction  is: 


—  5*  sin  a 
V*aVn  +  8*  cos  a 


(8) 


Summary  of  entire  calculated  results  leads  to  the  fol- 
owing  conclusions: 


1)  For  given  7. Me,  and  f»/to  ,  the  combination  parameter 

x 

DELTAPHT  varies  almost  linearly  with  ^  .  Therefore, 

the  least  square  principle  can  be  used  to  linearize 
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calculated  results:  DELTAPHT  -  K  ~y  where  K  is  slope 

2)  For  a  given  t  /t  and  r  as  a  parameter,  slope  K 
varies  with  M„  approximately  linearly,  i.e., 

XL 

K  (y,  M ,')  *•-  a  (y)  +  b(y)Ma 

where  a(y)  is  the  distance  between  the  origin  and  the 
intersection  of  the  line  and  axis,  *(y)  is  the  slope, 
and  a(y)  and  ,6(y)  are  all  functions  of  y 
In  addition,  a(y  )  and  b(  y)  also  vary  approximately 
linearly  with  y  ,  i.e., 

a  (y)  -  a0  +  ay,  b  (y)  «*  b0  +  bya 

3)  In  terms  of  engineering  applications,  in  order  to 
satisfy  the  viscosity  correction  for  a  nozzle,  the 
displacement  thickness  5*(a0  can  be  directly  added  to 
the  corresponding  coordinate  of  the  theoretical  nozzle 

wall,  i.e., 

*•(*)“  +$•(*). 

with  an  error  or  less  than  0.7%. 

4)  The  displacement  thickness  of  the  boundary  layer 

located  at  the  nozzle  throat,  8r,  ,  is  also  estimated. 

In  calculation  range,  8J  are  all  small, 

8;<0.0062A* 

In  practical  design,  the  effect  of  the  throat  boundary 
layer  may  not  be  considered. 

Based  on  the  foregoing  considerations,  tedious  calcu- 


lations  on  C-R  relation  including  pressure  variation, 
heat  exchange,  compressible  boundary  layer  of  laminar 
flow  can  be  substantially  simplified.  It  is  also 
possible  to  organize  the  entire  calculated  results 
into  approximate  formulas.  In  a  given  range  of 
parameters,  all  displacement  thicknesses  can  be 
calculated  with  those  formulas. 


The  approximate  formulas  for  the  combination  parameter  of 
displacement  thickness  after  arrangement  are: 


DELTAPHT= 


S*  /  Poh'tV*  \l/t 


-*■(%  Ma)^r 

K  (r,  Mg)  —  a  (y)  +  b  (y)  Mg 

•*  a0  4-  ai7  +  (bo+biy)  Mg 


(9) 

(10) 


For  a  laminar  current  heat-conduction  wall,  t  /t 

wo 

K  x  10* -=0.0088  +  0. 2ll4y 

+  (0.1131r-0.0041)  Mb 

For  a  laminar  insulated  wall,  t  /tQ  =  1, 


0.4, 

(11) 


K  X 10* -0.7173 +0.10647 

+  (0. 30327  -0.16G7)iW, 

The  error  introduced  by  the  approximate  formulas  is  less 
than  -  0.2%.  It  is  generally  less  than  -  0.003  mm  when  converted 

to  coordinates  of  the  nozzle  wall.  This  is  tolerable  in  engineering. 

The  appropriate  range  for  approximate  formulas  are: 
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y-1.33~1.4. 


For  heat  conduction  wall,  it  is  limited  to  tw/t  =  0.4.  Nozzle 
coordinates  are  calculated  by  the  Foelsch  method. 


REFERENCES 


[  1  ]  Gerry  E.  T.;  JE EE,  Spectrum.,  1970,  7,  51. 

[  2  ]  Anderson  J.  D.;  AIAA  Papor,  89-508. 

[  3  J  Anderson  J.  D.{  ACTA  Astronautical,  1975,  I, 
11. 

[  4  ]  Brunner  M.  J.j  AD-783218.  •  '  . 

[  5  ]  Anderson  J.  D.,  Harris  E.  L.;  AIAA  Paper  73- 
143. 

[3]  Greenberg  R.  A.  et  al.\  AIAA  J.,  1972,  1494. 

1 7  ]  Kuehn  D.  IS.;  Appl,  i’hyt.  Lett-.  1972,  <1,  3. 


[8]  Yan  Hai-shing,  Chen  Li-yin;  Chinese  Journal  of  Mechanics, 

1978,  No.  4,  274. 

[9]  FoeUch  K.;  North  American  Aviation  Report 
No.  NA-48-235,  Mar  1948. 

[10]  Crown  J.  0.;  NAOA,  TN-1651. 

[11]  MonslerM.  J.;  AIAA  Paper.  71-24." 

[13]  Stollery  J.  L.,  Smith  J.  E.;  J.  of  Fluid  ilttOu, 

'1982,  11,  225.  >'  . 

[13]  Simone  G.  A,;  AIAA  Paper,  72-709.  .  , 

[14]  Clark  P.  O.j  AIAA  Paior,  72-708. 

[15]  Kusnll  D.  A.jAIAA  Paper,  74-223. 

[19]  Wagner  J.  L.j  AD-781877. 

[17]  Coboo  C.  a,  Rvskotko  E.;  NACA  Rep  1294. 

A 


33 


